Exact Result on Hole Dynamics of 
One-Dimensional Supersymmetric t-J Model 
with Long Range Interaction 
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We obtain an exact result on single hole dynamics in one-dimensional Mott-insulators through 
the study on the supersymmetric t-J model with long range interaction at half-filling. We calculate 
the expression Gi s ih(x,t) for one spinon one holon contribution to the hole propagator. We find 
that Gi s ih(£ = 0,t = 0) = 3/8, which amounts to 75% of the spin-fixed mean particle density 1/2. 
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Single hole motion in Mott insulators is a nontrivial problem in strongly correlated electron systems. This issue 
becomes particularly intriguing in some one-dimensional systems: a spin liquid ground state causes a novel effect on 
hole motion; the underlying elementary excitations such as spinon (with spin 1/2 and no charge) and holon (charge 
e and no spin) provide a proper description to the hole dynamics. This situation seems to be realized in SrCu02 JO; 
in the angle-resolved photo-emission experiment on that material, the authors of ref. |l| identified dispersive peaks 
in the spectral intensity as holon and spinon dispersions. Meanwhile, the dynamics of the U — > oo Hubbard model 
has been studied intensively and important aspects of the spectral weight such as power-low singularity along the 
holon band || and the shadow band structure || have been clarified. However nontrivial is to what extent these 
properties survive in the t-J model with finite J. Further no simple analytical expression for the spectral weight has 
been obtained. Thus an analytical study on the hole dynamics in one-dimensional Mott insulators with finite J would 
be worthwhile. 

As a solvable t-J model, we consider the one-dimensional supersymmetric t-J model with long range interaction 
0. This model has simple but nontrivial physical properties; spinon and holon form an ideal gas |5j obeying the 
fractional exclusion statistics ||; at half-filling, this t- J model reduces to the Haldane-Shastry (HS) model • The 
dynamical spin correlation function of HS model has been derived in an extremely simple form Q . Thus this model 
is eminently suitable for our purpose. 

For this t- J model, Ha and Haldane |l(J studied the hole dynamics away from the half-filling and obtained the region 
of non- vanishing spectral weight of the spectral function in the frequency- momentum space ("compact support"). 
In addition, they identified excitation contents as simultaneous excitations of three spinons, three holons and one 
antiholon. At half-filling, however, the antiholons are suppressed. As a result, not only three spinon one holon (3slh) 
state but also one spinon and one holon (lslh) state contribute to the hole dynamics Eg] . 

In this paper, we consider the spectral weight of the hole propagator, which has been left as an unsolved problem 
in fiofl . To be precise, we calculate the exact expression for the (lslh) contribution Gi s ih(aj, t) to the hole propagator 
at half-filling. Our calculation has two steps; first we map the t-J model into a multicomponent Calogero-Sutherland 
model flll| , p"2| ■ Next we calculate the propagator in the resulting model with the use of the theory of the nonsymmetric 
Jack polynomials 113 14 and their variants 15 lq|. In the following part, we first present our result and then explain 



the derivation. 

We consider the t-J model 
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for the one-dimensional chain with even number of site N. Here V is the projection operator onto the Fock space 
without double occupation, ct (cj jCr ) and Si denote the electron creation (annihilation) operator with spin a(=], J.) 
and the spin operator at site i E [0, N — 1], respectively. The transfer energy tij and exchange energy between 



the sites i and j are related to each other by the supersymmetric condition: £y ; 
(N/n) sin [n(i-j)/N]. 

The hole propagator at half-filling is defined by 

G(x,t) = ((j>m\cl a (t) c o ,<t(O)|0hs)/(0hs,</>hs), 



Jij/2 = D i - . Here is given by 



(2) 



where |</>hs) denotes the singlet ground state at half filling, cj, a (t) is the Heisenberg representation of c) xa . The 
spatial separation x takes an integer value. We obtain the following expression for G\ s \\i{x, t) in the thermodynamic 
limit: 
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(-i) x r r~ q * In - q s 

G lslh (x,t) = / dq s / dgW -exp [-it (e h (g h ) + e s (g B ) - //)] cos [(g h + &) , (3) 

Jo Jo V ft 

where eh(<?) = (9 — 7*"/2) 2 and Csfa) = 9(1" — g) denote the dispersion of holon and spinon, respectively. The chemical 
potential ^ is given by — 7r 2 /12 [Q. The expression (g|) is our main result. 

The spectral function A{uj,Q) is defined by G(x,t) = (2n)- 2 j^duj dQe~ iut+iQx A(lj : Q) . In A(u,Q), the 
(lslh) contribution Ai 8 ih(o;, Q) for Q G [0, %] is given by 



[(a/ - Ch(Q)) (e B (Q) + tt 2 /4 - w')] / (Q + tt/2) 
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7rQ y w' — £h(<3) 

with u/ = a; + /i. Here 0[x] denotes the Heviside step function. The compact support of Ai s u 1 (o;, Q) is shown in Fig. 
1 by the dark region. In addition, the support of the full spectral function A(u>, Q) can be derived as the half-filling 
limit of the Ha-Haldane result JlCj]; the full support is also shown in Fig. 1 as the joint area of dark and light regions. 
The function Ai s ih(w, Q) has the inverse square root singularity along the curve lo — £h(Q) — fJ- indicated by the bold 
one. This singularity corresponds to the "main band" in the U — > 00 Hubbard model |^,|3|]. In our model, however, 
there is no equivalent of the "shadow band" within the one spinon one holon contribution. 

Next we consider the static limit. For Q S [0, 7r], the momentum distribution function is given by 

/oo 
dxe~ iQx G lslh (x,t = 0) = it' 1 AxcSmy/l - Q/tt + tT 2 VQ ¥ - Q)- 
-00 

In p!^| , on the other hand, the full result was found to be G(Q, t = 0) = 1/2. In Fig. 2, we show Gi s ih(Q, t — 0) and 
G(Q,t — 0) = 1/2. We can see that the (lslh) contribution is dominant near Q — while most contribution near 
Q = 7r comes from (3slh) states. Additionally, we find that Gi s ih(^ = 0,t = 0) = 3/8. This means that the (lslh) 
contribution amounts to 75% of the mean density 1/2 with fixed spin. 

Now we present the derivation of (^). We denote a state vector \<f>) without double occupation by 

10)= tfM.M) n 5 r n h \\ p )- w 

{x},{y} ie{x} je{y} 

Here {x} = {xi, • • • , xm} denotes the set of coordinates for M down-spin electrons and {y} = {yi, • ■ • , yp} denotes 
those of P holes. The state \F) is the fully-polarized up-spin state. The operators — c\ , Cjj and h\ = Cjj represent 
the spin- lowering and hole creation operator of site i, respectively. 

The wavefunction ip ({x} , {y}) is defined only when each argument takes an integer value. However we can regard 
it as a restriction of a complex function $ ({X} , {Y}) of {X} = {Xi, ■■■ , X M } € C M and {¥} = {Yi, ■ ■ ■ , Y P } £ C p 
if the relation ({x} , {y}) = $ ({X} , {Y}) hold for X k = exp(\2irx k /N) with k E [1,M] and Yi = exp(i2Tryi/N) with 
/ G [1,P]- We can impose further two conditions on $ ({X} , {Y}) without loss of generality: 1) symmetric under 
pairwise interchange of X. 2) antisymmetric under pairwise interchange of Y. 3) zero when X k = Xi, Y k = Yi or 
X k =Y h 

From now on we concentrate on the SU(2) Yangian highest weight states (YHWS) (l^Jl^l because the (lslh) states 
belong to those states. The YHWS are the states with polynomial-type wavefunctions where the degree of each 
variable is in the range (0, N). The YHWS form a subset of the Fock space and the residual states can be generated 
by successive actions of a Yangian generator on YHWS [pi. For YHWS, it is known ]2(J that wavefunctions are given 
by those of the SU(1,1) Calogero-Sutherland model 



• u + /> ' d A \ 2 x . 1/.,/. 



i=l v ' i<j 

e$. (5) 



Here we denote Z = (X, Y) = (Xi ■ ■ ■ Xm^Yi • ■ ■ , Yp)- Sij is the coordinate exchange operator of the pair i and j. 
The eigenvalue e in (|[) is related to E in ([j]) as 

E = 2n 2 (e + SM:P)/N 2 , (6) 

with e M ,p = N(N 2 - 1) /12 - (Y 2 + 2)P/4 - Y 2 Af/4. Wavefunctions for YHWS have the form 
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M M+P M+P 

£<j i<j 2—1 



(7) 



Here &(X, Y") is a polynomial symmetric with respect to the interchange of (Xi,Xj) and that of (Y, Yj). The degree 
of Xi in $ is less than or equal to Y sp = N — 2M — P ^0|. We can regard N sp as the number of spinon excitations. 
The relation 5 t z ot = N sp /2 holds only in YHWS. 

Now suppose that a hole with S z = 1/2 is injected to the singlet ground state at half-filling. Then the states with 
S% oi . = 1/2 and P = 1 are relevant intermediate states; P = 1 means that one holon is excited. In the intermediate 
states, only (lslh) states (N sp — 1) are YHWS because 2S't ot = 1. In the following, we diagonalize the Hamiltonian 
(§) for M = M = N/2- land P = 1. 

First, we define the basis function ® JK for J e [0,M ] and X e [0,M ] as *, 7A ' = F J e_ftr(X)$ - Here e K (X) 
is the elementary symmetric function of order K, which is defined by J2i<k 1 <k 2 ---k K <m„ IT;e{fc} X% with {k} = 
{fei, k2, • • • , fcif }■ The function <I>o is defined by 

M M n M 

$ =yn^n(^- x i) 2 n^- y )' ( § ) 

i— 1 i<j i— 1 

Now we consider the case J + K < Mo- In the present basis, the diagonal element j k\H\$> j k ) is given by 
ejk = 2 J 2 - 2JM - 2Jr 2 + 2KM + J + K - M + e H s with e H s = (Y 2 /4 - l)N/6. The off-diagonal element 
(^j+ijif-ilWl^^if) is (~ l) L 2(Mo — J — AT) for L > and zero for L < 0. The Hamiltonian matrix in this basis is 
triangular and hence it can be solved recursively. As a consequence, we obtain 

K 

®JK = ^ a L^J+L,K-L, (9) 
L=0 

for J + K < M . Here the coefficients a L are given by a = 1 and a L = (-l) L+1 2" 1 7r" 1 / 2 r [L - 1/2] fL\ for L > I. 
For J + K > Mo, the eigenfunction is obtained similarly as 

Mo— K 

*JA' = XI a L^J-L,K+L- ( 10 ) 
L=0 

with the eigenvalue e JK = 2 J 2 - 2JM Q - 2K 2 + 2KM + M - J - K + e H s- 

With these preliminaries, we calculate the (lslh) contribution to the hole propagator 

Gi s ih(a;,t) 

M M 

= J2J2 I<^|co,iI^hs)| 2 /((*jaH*jk)(*hs|*hs)) 

J=0K=0 

(-If exp [-it (E JK - Ens) +i2n{J + K+ 1) x/N] . 
Here the ground state wavefunction $hs at half-filling is given by 0J^] 

N/2 N/2 

$n S = l[(X i -X j ) 2 l[X i . (11) 

i<j i=l 

The eigenvalues E JK and i? H s of \®jk) an d |$hs) for Ha i n (@) are given by 2n 2 (ejK + £m .i) /N 2 and 
2-7T 2 (ehs + £m +i,o) /N 2 , respectively. In the following, we consider the norm first and then discuss the matrix 
element. 

For general M and P, the inner product between |$) and |$'} can be written as 

*') = N M+P $']„ /(M!P!), (12) 

when the wavefunctions &(Z), Q'(Z) are polynomials of Z. Here we define a scalar product [$,<!>'] on polynomials 
by the constant term in $(1/Z)$'(Z); the symbol 1/Z denotes {1/Z\, • • • , l/Z M +p)- 
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The norm [$hs, < J ) hs]o for ( |Tl] ) has been known as A N — 2 N / 2 N\ [El 22|. In addition, the norm [$jk, ®jk]o for 
@ or (|l~o| ) is obtained as 

Aat r [JV/2 - X - A JA -] r [K + 1/2 + A JK ] 

JV r[iV/2-ir+l/2-A JK ]r[K + l + A JK ]' 1 j 



from the theory of a variant of non-symmetric Jack polynomials |16|. In Jig), Akj is 1/2 for J + K > Mo or zero 
otherwise. 

Next we consider the matrix element (&jk\co,i\&hs) = ( < f ) Ji<"|^o'5'o"l^ ) Hs)- We present the derivation only for 
J + K < Mq because the following discussions are applicable also to the case J + K > Mo . 

Now we introduce homogeneous symmetric polynomials J^{X) = J filjtl2 (X) of X = {Xi, ■ ■ ■ , Xm„} with /i = 
G"i,M2) by 

Mo K 

e K (X)e K ,(X)= E 

8k+k' ,iii+ii 3 du Jn(X) (14) 

fii—K iii— 

for K,K' £ [0,M ]. Here <5 aj 6 denotes the Kronecker's delta and the coefficient is given by 

r [mi - X + 1/2] T[K- /X2 + 1/2] r [ M i - + 1] 

r [1/2] r [mi - /v + 1] r [x - M2 + 1] r [mi - M2 + 1/2] ' 

The polynomials J M are called the symmetric Jack polynomials [ p3| . They are useful in our calculation because they 
form an orthogonal basis with the following norm |2^] : 

[J,*HS, J,$HS]0 - ^727 F 2 [M1-M2 + 1] 

A T[pLi + (3- i) /2]F [Mq -m + i/2] 

f = \T[ f x i +2-i/2]T[M - f i i + (i + l)/2Y [ ' 

Here $ HS denotes nf<° (** - Ilf=i 

The wavefunction of /io| < I ) jjf) is given by ^Jir({-Xi}j = °i , Y = 1). With the use of (|l4|). this wavefunction can be 
written as 

if M K-L 

E E E(- 1 ) i ^ 1 ^ 2 ^<" Lj ^Hs, (i6) 

up to a phase factor. 

On the other hand, the wavefunction of Sq |3>hs) is given by j|] nf=i (^» — I) 2 ^hs- Further the decomposition of 
this wavefunction in terms of J m 4>hs is known as pl| , p2] | 

E i-ir^B^j^^s, (17) 

with B^.^ = T [1/2] T [mi - M2 + 1] /r [mi - M2 + 1/2]. From the expressions @, |[|), @ and (0), we obtain a 
simple expression for the matrix element 

\(^j K \hlS+\^ s )\=g N T[K + l/2]/T[K+l] (18) 

with g N = tt- 1 I' 2 {N/2) n ' 2 ' 1 {N - l)!/(iV/2 - 1)!. 
A similar calculation for J + K > Mq gives 

I($,/a-|^5*o + I $ hs)I - 3Arr[Af - K - l/2]/T[M - (19) 
As a result of (§J|), © and ©, we obtain the expression for the hole propagator 
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G lalh (x,t) = 1/(2N) 

4 (-1)" N/2 ^- 2 r [AT/2 + 1/2] r [g + 1 /2] 

+ irN 2 r [N/2 - K) T [K + 1] 

e - WB ^cos(i27rx(J + A' +1)/JV), (20) 

with £' 7X = 2tt 2 (2 J 2 - 2 JM - 2i^ 2 + 2KM + J + K - M - 1/6) /TV 2 + tt 2 /3. Further wc introduce q h = 2ttJ/N 
and q s = 2nK/N and take the thermodynamic limit N — > oo. Finally, the expression ( p0| ) reduces to (^). 
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FIG. 1. Region of non-vanishing spectral weight function of the hole propagator for Q £ [0, n]. The dark region represents 
the contribution from states with one spinon and one holon. At the lower edge indicated by the bold line, the spectral weight 
has the inverse square root singularity. 

FIG. 2. The contribution from states with one spinon and one holon to the momentum distribution function Gi s ih(Q,i = 0) 
for Q e [0, tt]. Full result G(Q, t = 0) = 1/2 is also shown. 
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